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Abstract 

Collapse of the rotating spheroid is approximated by a system of or- 
dinary differential equations describing its dynamics. The gravitational 
^ ' potential is approximated by the one of the iniform Maclaurin spheroid. 

OO , Developement of gravitational instability and collapse in the dark matter 

medium do not lead to any shock formation or radiation, but is character- 
ized by non-coUisional relaxation, which is accompanied by the mass and 
angular momentum losses. Phenomenological account of these processes is 
done in this model. Formation of the equilibrium configuration dynamics 
I of collapse is investigated for several parameters, characterizing the con- 

■~«.^ ■ figuration. A very long gravitational wave emission during the collapse is 

i-C ' estimated, and their possible connection with the observed gravitational 

lenses is discussed. 

o 

H 

to '. 1 Introduction 

a 

A study of of a formation of dark matter objects in the universe is based on 
N-body simulations, which are very time consuming. In this situation a sim- 
' plified approach may become useful, because it permits to investigate rapidly 

. by analytical or by simple numerical calculations many different variants , and 

obtain crudly some new principally important features of the problem which 
could be lost or not visible during a long numerical work. 

The modern theory of a large scale structure is based on the ideas of Zel- 
dovich (1970) about a formation of strongly non-spherical structures during 
non-linear stages of a development of the gravitational instability, known as 
"Zeldovich's pancakes". The numerical simulations started by Doroshkevich et 
al. (1980) had been performed subsequently by many groups, revealing compli- 
cated structures of the new born objects (see e.g. Doroshkevich et al., 1999). 

Here we derive and solve equations for a simple model of a dynamical be- 
haviour of a compressible rotating spheroid in which a motion along both axies 
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takes place in the common gravitational field of a uniform spheroid, and un- 
der the action of the isotropic pressure (taken into account in the approximate 
non-differential way), and otherwise independently on each other. Due to non- 
interacting particles of the dark matter a strong compression docs not lead to 
a formation of the shock wave and additional energy losses due to radiation at 
increasing of a temperature.. All losses are connected with a runaway particles, 
and relaxation consists only in the phase mixing leading effectively to a trans- 
formation of the kinetic energy of the ordered motion (collapse) into the kinetic 
energy of the chaotic motion, and to creation an effective pressure and thermal 
energy. In presence of the chaotic motion the pancake formed during the col- 
lapse has a finite minimal thickness, because particles cross the equatorial plane 
non-simultaneously, forming therefore the effective pressure, which "stops" the 
contraction. These effects may be described approximately in the hydrodynam- 
ical approach, in which there are no shock waves and the main transport process 
is an effective bulk viscosity, leading to relaxation and to a damping of the or- 
dered motion. This relaxation is based on the idea of a "violent relaxation" of 
Lynden-Bell (1967). The system of the ordinary differential equations is derived, 
which describes the dynamical behaviour of the compressible spheroid, where 
the relaxation and losses of the energy, mass, and angular momcntiim are taken 
into accound phenomenologically. In absence of any dissipation these equations 
describe a self-consistent conservative system, where gravitational and thermal 
pressure forces are present. Namely, the equilibrium solution of these equa- 
tion describes a Maclaurin spheriod, where the density is not prescribed but is 
found self-consistently, when the effective entropy is known. We consider here 
only non-relativistic dark matter with a relations between the density p, energy 
density E, pressure P, and specific entropy S as 

Q 

pE=-P, P ~ p^/^, E ~ p^/^ at constant S. 

Solution of simple ordinary equations reveals some interesting features of the 
collapse of non-interacting matter, which had not been noticed earlier. The 
collapse of a sphcrioid, leading to the formation of a pancake may be continued 
by the fornation of a transient oblate figure, strongly elongated along the axis 
of a symmetry with a ratio of the axis c to the equatorial radius a : c/a ^ 5. 
Formation of such transient figures, appearing only at low rotation, could be 
expected in central parts of the collapsing dark matter objects during initial 
stages of the collapse, when the influence of the outer non-uniform part of the 
cloud is still negligible. 

We have found, that at a realistic rate of the relaxation rate with a character- 
istic time equal to 3 local Jeans times, the oscillations do not damp completely, 
and about 1/10 of the initial amplitude of a velocity survives after at least 10 
oscillationl periods. That means that the most massive dark matter objects 
may be still in the oscillatory state. These oscillations, as well as variable grav- 
itational fields of the collapsing dark matter objects could be important for the 
the interpretation of the observed picture of the cosmic microvave background 
(CMB) fiuctuations. Very long gravitational waves emitted mainly during the 
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first stage of the pancake formation (see Thuan and Ostriker, 1974; Novikov, 

1976) could be also important for the CMB fluctuations (Doroshkevich et al., 

1977) , as well as for the gravitational lensing of the distant objects. 



2 Properties of a Maclaurin spheroid 



Let us consider a spheroid with semi-axics a — b > c 

^+4=1, (1) 

and rotating uniformly with an angular velocity n. Let us approximate the 
density of matter p in the spheroid as uniform. The mass m and total angular 
momentum M of the spheroid arc connected with the angular velosity and semi- 
axies as (Landau, Lifshits, 1988) 
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Introducing k — c/a < 1, we may express the gravitational energy Ug as 

ZGm? arccosfc 

Un — 



(2) 



(3) 



5a Vr^rfc2' 

what gives U — —iGm^/ba for a sphere. The gravitational forces inside the 
uniform spheriod Fx-, Fy, Fz are expressed via the grvitational potential </> as 
F = — V0, and are defined by the relations 



2T:Gpx 
l-fc2 



k arccos k 



(4) 
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k arccos k 



The gravitational potential with normalization </> - 
spheroid is expressed as 

(j) = -i(xF^ + yFy +zFz)- 2TTGpa 



-> at infinity, inside the 
, k arccos k 



During the process of relaxation a form of the spheroid changes between 
oblate and prolate ones. When c > a = b the following relations are valid 
instead of (|)-(|) 
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(7) 
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^ Vk"^ — I 

Here the positive branch of Archfc should be used. 



-Tz = -m — 



3 Equations of motion 



For a real Maclaurin equilibrium noncompressible spheroid the oblateness k is 
connected with the angular velocity (or angular momentum) 



fc(l + 2A:2) 



arccos k 



3fc2 
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1 - fc2 27rG/9 
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(9) 



(1- fc2)3/2 

We consider a case of non-collisional particles when the density varies in time, 
and approximate the form and gravitational potential of the body by the form 
and gravitational potential of Maclaurin spheroid. Here pressure does not ne- 
sessary balance the garvitational and centrifugal forces. Consider a uniform 
spheroid with a constant mass and angular momentum, and with a total "ther- 
mal" energy of non-relativistic dark matter particles Eth ^ F~2/3 ^ (a6c)~2/3 
{V is a volume of a spheroid, a = b). Equations of motion for this spheroid are 
written for the radial accelerations at the equator and at the pole for semiaxes 
a, c 



25 M2 3 Gm 

Tm^a^ ^ 2a2(l-fc2) 

Gm ^ 
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(11) 



Here Eth,in, o-in ~ bin, Cin are the initial values of corresponding parameters. 
Equations (|o|),(|ll|) describe the dynamics of the concervative system with a 
linear dependence of the velocity on the coordinates 



ax by cz 

abc 



and Lagrange function 
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rot -^Ug + Eth, 
^ — Eth^ini^inbin^in^ ^ 



(12) 
(13) 



(a6c)2/3 ■ 

The function Ug defined in (||), and the rotational energy from (|l|) should be 
used during the variation of L in the form 
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5(a + 5) VT^fc^' 
The kinetic energy is equal to 



k' = 
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(14) 




dx dy dz 



m 1-2 



(15) 



Equations of motion ([l0|) , (ITTl) are the Lagrange equations with the Lagrange 
function ([r^)p] Solutions of these equations describe either pure oscillations, 
or the total disruption, depending on the initial conditions. For the prolate 
spheroid the following relations are valid instead of (p^ , (|Tl|) , (p^ 
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c = 3 
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(16) 
(17) 
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4 Equations of motion with dissipation 



(18) 



In the reality there is a relaxation in the coUisionless system, connected with 
a phase mixing which called "violent relaxation" (Lynden-Bell, 1967). This 
relaxation leads to a dissipation of the energy of the kinetic motion and increase 
of the chaotic (thermal) energy and pressure. As a result of this dissipation the 
kinetic motion will suffer from an effective drag force, which account is described 
phenomenologically by adding of the terms 



and 



'^rel 



'^rel 



(19) 



^We consider here for simplicity that the relaxation is accompanied by the isotropization 
of the distribution function in the velociyu space. The violent relaxation has a non-collisional 
origin, and if the anizotropic instability has a lower increment the anizotropy in the velocity 
space may be preserved. In the limiting case of a pure anizotropic relaxation we have three 
(or two in our case) pressure and entropy functions, so that 
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2 £z 



We should have the terms lOSx/a^m and lOfz/c'^m instead of the last terms in (px|),(pJ^, 
respectively. A corresponding changes should be done also in the equations determining 
increase of entropies, and different losses. 
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in the right parts of equations (]lO|),(pl|) respectively. We consider the nonrela- 
tivistic dark matter with the relation between the pressure P and the thermal 
energy Eth of the whole spheroid as 

E.. = l^. (20) 

Dissipation (|l^) leads to a heat production. Let us find a rate of a heat pro- 
duction leading to the growth of the entropy of the spheroid matter. Write 
equations of motion (|lo|),(|ll|) with account of (|5|),(13),([l9|) separately for a, b 



and c, allowing also the energy, mass and angular momentum losses 

.dUg .dUrot .dEth 1 a' 

---a— a— a— -m , (21) 

oa oa oa 5 Trei 

■dUg -dUrot ;dEth 1 

Ob Ob Ob 5 Trei 

.dUg .dUrot .dEth 1 c2 
= -c— ^-c— c— -m , (23) 

OC ac OC 5 Trel 




Trel = 2'Karel\j (^4) 

Note that in (pT|)-(p^) the gravitational energy should be taken in the form (p^. 
We have scaled here the relaxation time r^eZ by the Jeans characteristic time 



2-K 2-K abc , , 



with a constant value of arei- Collecting together (|2l|)-(23) we find the energy 
balance relation in the form 



dUtot _ dUg dm ^ dUrot dm ^ dUrot dM ^ dEth d£ ^ Ukin Ukm dm 
dt dm dt dm dt dM dt d£ dt r^ei rn dt ' 

Utot — Ukin +Ug + Urot + Eth- 

The process of relaxation is accompanied also by the energy, mass, and angular 
momentum losses from the system. We suggest, that these losses take place 
only in non-stationary phases, so these rates are proportional to the kinetic 
energy Ukin, the characterictic times for mass, angular momentum and energy 
losses Tml , tmi ,Tei are considerably greater then r^ei . The equations describing 
different losses may be phenomenologically written as 

dUtot Ukin 



dt Tml 

dm mUk' 



kin 



dt UnT, 



gTml 



Ug < 0, (28) 
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(29) 

dt UgTMl ■ ^ ^ 

Scaling all the characteristic times by the Jeans value, like in the case ( p^ of 
Treh we have 



Tel — CtelTj, Tml = CXmlTj , TMl = OLMlTj, (30) 

with constant values of (i = el, ml^ Ml). Combining (p^) , (p9|) , we find a 
relation 



M 



m \ "Ml 



(31) 



where TOi„ and Mi„ are the initial values of corresponding parameters. The 
function £ determines the entropy of the matter in the spheroid, so that, using 
Isl), we have 



dEth^pdV 
dt 



1 d£ 



th 



dt [abcy/^ dt 
, V 



(32) 
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Using (p^-(p9|) in (p6|) we obtain the equation for the entropy function £ in 



presence of different losses in the form 
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5 Non-dimensional equations 

The equations, describing approximately the dynamics of the fomation of a sta- 
tionary dark matter object include equations of motion ( p^, ( [ll]) with adding 
terms from jl^); energy equation (^) with account of (|^),(|l3|)7p5|); and equa- 
tions (27)-(P9^ describing the losses of the energy, mass and angular momen- 
tum. The mass losses are also taken into account in the equations of motion, 
written in the form (pll)-(p3|). For obtaining a numerical solution of these equa- 
tions let us write them in non-dimensional variables. Introduce the following 
non-dimensional variables 



m 
mo' 



M 



(34) 
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- Il 

Po Ua Ua to to 

The process, with account of mass and angular momentum losses, is described 
by the following non-dimensional system of equations 



~ t a c 

t = — , a = — , c = — , ? 
to ao ao 
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Ua Ua £o 
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a a{2a^ + c^) 25 3 to 

° ^ S ^ 1 2~3 + o~°^lfc 
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for the oblate spheriod a = b > c, and 



a*; = 
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k- ak 



a2t = 



1 — fc Ofc 
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(36) 
(37) 



(38) 
(39) 

(40) 

(41) 

k a2k. 



for the prolate spherfod with a — b < c. In both cases + alk 
The rotation phase / is calculated in ( p^ ) for comparison of the rotation and 
oscillation processes. When the spheroid is close to the sphere |fc — 1| <C 1 we 
have the following expansions 



~ 4 1 . 

7 k^ 
a-k ~ — ;ri alk 
6 6 
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(42) 



1 



^, a2k 
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— at a = 6 < 



fc+1 6 k+1 6 

In the case of a very compressed spheroid (pancake) we have the following 
expansion 



flfe = fc H — k k^, oil, 

2 4 3 ' 



-!!: + 2fc- — fc^, a2fc = 1- -fc + 2fc2. (43) 



In equations (|35|)-(|39|) only non-dimentional variables are used, and "tilde" sign 
is omitted for simplicity. The non-dimensional relaxation times in these 
equations are written with account of (E5|),(PG) as 



27rQ!,- 



(44) 



8 



The scaling parameters 

to, ao, mo, Mo, poi Uq, flo, So 
are connected by the foUowing relations 



— ^, M2 = Gaom^, C/q = ^, po = ^fo = 



(45) 

so that in non-dimensional variables m — Ana^cp/i, M = 2ma^il/5. The total 
nondimcnsional energy Utot — Utot/Uo is written as ("tilde" is omitted) 

m(2h^ + c^) 3m2 5 E 
10 5a 4ma2 (0^0)2/3 

6 Equilibrium configuration and linear oscilla- 
tions 

In equilibrium there is no dissipation and the following relations follow from 
(^5|) , (|3^) , connecting the equilibrium parameters 

25 3 m 10 £ 

4 m2a3 2 a2 " + 3am (a2c)2/3 ~ ' ^ 

At given m, M and £ these equations determine the configuration of the sphe- 
roid. Finding a from ( ^ ) 

^ (49) 



9m2 A:5/3 ' 



we obtain after its substitution into ( [47|) the relation for determining fc in the 
form 

„^+2.a2.-^i^^^!^!^;^Ht = o. (60) 

4 m 5 ^ ' 

Taking into account (^0|), we come to the relation 

A;(l + 2fc2) , 3fc2 75Af2 .,3 ^ 25 A/^ fc 

— ^arccosfc- = 1 = — — -. 51 

(1 — fe^)3/^ 1 — fe^ 4 mc 6 m'^ a 



The last relation in ( |5l| ) follows from the connection (U9|) between a and £, and 
in this form (|^) coinsids exactly with the equation (^, determining the form 
of the Maclaurin spheroid with 
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of a given mass, density and angular momentum. In our case the density is 
determined by the entropy £. For the spheroid with unit non-dimensional mass 
m and semiaxis a we obtain in equilibrium relations determining M and £ as 
functions of k 

£ = Afc5/3^2fe, ^-^{alk +2ka2u). (52) 

10 25 

The spheriod, described by ordinary differential equations for semiaxies a and 

c has two levels of freedom in oscilations corresponding to the first p- mode, 

and the first fundamental /-mode (see Cox, 1980). Consider (without loosing 

a generality) oscillations of a spheroid with m = a — 1. The £, Utot,n^ and M 

losses are quadratic to perturbations, so these values are conserved during linear 

oscillations. For the perturbations of da and 5c we have the following equations 

from M),m 



T5 3 £ 20 £ Sqj 

Sa = -—M^Sa + -S{alk) -3alk Sa ~ lOj^Sa - -^^5k - — , (53) 

£ ^Vi £ Sc 

Sc = -3(5(a2fe) + 6 a2fe Sa - lO^Ja - -^^5k - —. (54) 

The perturbations are expressed via perturbations Sa, Sc as 

S{ak) = -a2kSk, S{alk) ^ ^ ^ a2k+2kalk ^^^^ 

S{a2k) ? , 9 ^ Sk, Sk^Sc~Sa; Trei ^ 2Trarei\ t:- 
1 — k^ V 3 

The equations (^) , ([S^) with account of (^5|) may be solved analitically at weak 
damping for any k. Consider here for simplicity two limiting cases k = 1 
(sphere), and fc <C 1 (pancake). The sphere is related to the nonrotating body 
with M = 0. In this case we have also 

2 13 

Ofe = 1, alk = -g, a2k = -, ^ ^ J^^ (56) 

Siak) = — 5", S(alk)^—, d(a2fe = — — . 

3 5 15 



With account of (56) we have from (B3h,(p4) the following equations 



c" 1 r, c ^o, ■■ 13 ^ Sc 

Sa = Sk — Sa , dc = Sk — da . (57) 

15 '^rel 15 '^rel 

These equations are completed by relation Sk — Sc ~ Sa. Taking Sa, Sc ~ 
exp(— icji) we come to the characteristic equation 
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The solution of a characteristic equation may be obtained analytically for a 
weak damping at a; 3> 1 /TVeZ • The 4 roots of this equation are 

i [a i 

Wi,2 = ±1 - , W3,4 = ±W - ~ (59) 

ZTrel V ZTrel 

The first root corresponds to a pure radial oscillations, and the second one 
describes the fundamental mode which does not change the volume and which 
is preserved in the non-compressible fluid. In fact, for a pure radial oscillations 
with ^fc = we have Sc = Sa, and the equations ( |57| ) reduce to Sa — —Sa — 
with a characteristic equation + — 1 = 0. The solution of this equation 



is wi.2 = ~2t~i -'- Y -'^ ^ 47^' "^tiich for a small damping obviously correspond 
to the first two roots in (Isof). The damping increments are equal for these two 



modes, what may be seen also from the numerical solution of (|35|)-(|39[) for a 
general nonstationary problem at large t (see next section). Note, that in the 
radial oscillation mode Sa = Sc, and in the fundamental mode Sa w —O.bSc. 
For the case of a pancake fc <C 1 we have from (E3),(p2h, (pq) 



a, = |, aU = a2, = 1, £ = Afc5/3^ ^,^2 ^ |^ (g^) 

S{ak)^-Sk, S{alk) = 2Sk, S{a2k) = -^. 



With account of (60) we have from (b3[),(p4) the following equations 



Sa — Sa — 3Sk , Sc — — — Sk — 3Sa , (61) 

4 Trel k Trel 

with Sk = Sc — kSa. Taking Sa, Sc ~ exp(— iwt) we come to the characteristic 
equation 



The solution of a characteristic equation may be obtained analytically for a 
weak damping at w 3> 1 /Trei ■ The 4 roots of this equation are 



^^1,2 = 




In the pancake the oscillation modes (1,2) and (3,4) are almost independent 
on each other. In the first two modes Sa ~ kSc <^ Sc, what corresponds to 
rapid oscillations of the thickness of the pancake at almost fixed radius. Two 
other modes with a low frequency are characterized by relation Sc ^ kSa <^ Sa, 
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and correspond to the oscillation of the equatorial radius of the pancake at 
approximately the same thickness. Here both modes have the same damping 
increments because we have chosen the same relaxation time for the damping 
of two degrees of freedom, which may differ in reality for the flat pancake. 

The model, developed here may be easily generalized onto the 3 dimensional 
ellipsoids, where instead of analytical formulae for the gravitational potential 
and forces we'll have elliptical integrals. This model will have three degrees 
of freedom, with the additional one, corresponding to non-axisymmetric per- 
turbations. Contrary to the axisymmetric modes, which are always stable, the 
non-axisymmetric mode of the spheroid will become unstable at sufficiently large 
M, first due to a secular instability appearing only at presence of damping, and 
at larger M the spheriod will become dynamically unstable independently on 
the presence of the damping. In that case we may expect that in presence of 
damping the dynamical instability will appear at larger M, than in the pure 
conservative case (see Chandrasekhar, 1969). We suppose to investigate the 
three-dimensional model and its stability in the near future. 



7 Numerical results 

The development of gravitational instability starts from almost Hubble expan- 
sion velocity which is transformed into contraction as a result of a developement 
of gravitational instability. To simplify the problem we start from a nonlinear 
phase of the instability with an expansion velocity smaller than the Hubble 
value. We start the simulation from the spherical body of the unit mass, small 
or zero entropy £in <C 1, and negative total energy Utot < 0. We also specify 
the parameters, characterizing different dissipations Q!rei, cteZ, ctmZ, ciMi- The 
following parameters have been used in all variants of calculations 

flm = 1, Ci„ = 1, rain = 1, Otrel = 3, Uel = a^i = CtMl = 15. (64) 

Other four initial parameters di„, Ci„, M and £ used in calculations are listed 
in the Table 1, where also the initial Utot is represented, which is calculated 
from (H). 

We have used in the calculations the Runge-Kutta code from (Press et al., 
1992). The results of the calculations are very transparent, and are represented 
in figures 1-16. The odd number figures represent time dependences of a and 
c axies on time, together with the phase = sin(p/i) — sin(Jp fldt), permitting 
the comparison between the oscillational and rotational motion. In the even 
number figures the entropy £, current mass m and current total energy Utot are 
given as functions of the time for the choosen set of the relaxation parameters 
(|64[). Here the entropy behaviour is the most interesting, because it is based 
on the realistic estimation of r^ei ■ The characteristic times of other losses are 
more difficult to estimate, so they had been choosen equal to ST^ei, making these 
losses practically unimportant. 

After the initial expansion the collapse happens, during which the linear 
size decreases 10-40 times, and in the process of the subsequent relaxation the 
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N 




Cm 


M 


£ 


Utot 


1 


0.5 


0.5 


0.3 


0.01 


-0.4025 


2 


0.5 


0.5 


0.5 


0.01 


-0.2025 


3 


0.2 


0.2 


0.5 


0.01 


-0.2655 


4 


0.5 


0.5 


0.1 





-0.5125 


5 


0.9 


0.9 


0.1 





-0.3445 


6 


0.9 


0.7 


0.1 





-0.3765 


7 


0.99 


0.94 


0.01 





-0.3155 


8 


1.3 


1.1 


0.01 





-0.1409 



Table 1: Initial conditions and parameters taken in numerical solution of equa- 
tions (p5[)-(p9[), in addition to the fixed values from (64); and initial value of the 
total energy, according to (|2F 



object approaches the equilibrium, described by relations (p2[). The relaxation 
and entropy production, mass and energy losses happens in the most rapid 
stages at maximal compression, what is reflected in the steps of the functions 
£{t), m{t), Utot{t) , well visible in figures 8,10,12,14,16. The following remark- 
able properties should be noted. 

1. For a low initial angular momentum, which is probably characteristic for 
the collapse of the primodial dark matter clouds, the collapse leads initially to 
the "pancake", in accordance of Zeldovich (1970) theory. On the subsequent 
stage of the expansion the spheroid becomes prolate, and during subsequent 
phases the spheroid is oscillating on the fundamental mode, changing the form 
between the prolate and oblate ones. It may be seen especially well in figures 
11 and 15, where temporally the prolate figure is formed with c/a ^ 2 — 5. We 
may expect the formation of such transient prolate figures in the central parts 
of the collapsing dark matter cloud, where the density may be taken as almost 
uniform, and where the characteristic period of oscillations is much less than 
the time of the whole body collapse. 

2. The phenomenological relaxation time accepted in the calculations seems 
to be rather reahstic (Lynden-Bell, 1967), so the calculated process of the damp- 
ing should be close to a reality. We have obtained, that oscillations preserve 
during more than 10 initial characteristic (Jeans) times tj. That means, that 
in the most extended dark matter objects with the largest tj the damping of 
oscillations is still in progress, and the dark matter objects may be presently in 
a nonstationary state. The collapsing dark matter ofjects and their oscillations 
could influence the visible picture of the fluctuations of the cosmic background 
microwave radiation, which are now under intensive investigation (see Naselsky, 
Novikov, Novikov, 2002). 

3. The description of the dynamics of the spheroid by the system of equa- 
tions with two degrees of freedom give a possibility of the exitation of only two 
oscillation modes. In the case of a low rotation when two degrees of freedom 
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are strongly coupled the two modes reduce to the first radial p-modc, like the 
main mode of the radial oscillations of the sphere, and the fundamental /-mode, 
which almost does not change the total volume, and which survives in the non- 
compressible sphere (Cox, 1980). The lengths of two axies oscillate with the 
same period. Excitation of both modes may be distinctly seen in calculations. 
For some initial conditions the main remaining mode is the p-mode (figures 
7,15); for others (figures 9,11,13) the fundamental /-mode was mainly excited. 

For rapidly rotating bodies, with a large equilibrium ratio of a/c, the oscilla- 
tions in two axies are almost independent, and their lengths are oscillating with 
substantially different periods (figures 1,3,5). 



8 Emission of very long graviattional waves 

Gravitational radiation is produced during the collapse of the nonspherical body. 
Gravitational radiation during a formation of a pancake was estimated by Thuan 
and Ostriker (1974), and was improved by Novikov (1975), who took into ac- 
count the most important stage of the radiation during a bounce. The formula 
for the estimation of the total energy emitted during the collapse and bounce is 

Ugw « 0.01 ( ^] { ^ Mc\ (65) 



Here Oe^ is the equilibrium equatorial radius of the pancake, and Cmin is its 
minimal half-thickness during the bounce, Tg = 2GM/c^ is the Schwarzschild 
gravitational radius of the body. From our calculations wc have aeq/cmin < 100. 
The value of Vg/aeq we estimate using the observed average velocity of a galaxy 
in the rich cluster Vg 3000 km/s, and taking 

f!^)'«10-4. (66) 

aeq y C / 

Than the energy carried away by the gravitational wave (GW) may be estimated 
as Ugw ~ 10~^^Mc^. If all dark matter had passed through the stage of a 
pancake formation, than very long GW with a wavelength of the order of the 
size of the galactic cluster have an average energy density in the universe 

ecw « lO-^VdmC^ » 3 • 10-23erg/cm^ (67) 

Here we have used for estimation the average dark matter density pdm = 3- 10~^° 
g/cm^. The average strength Eqw of the very long GW may be estimated, 
taking the relation (Landau and Lifshits, 1993) 

where h is metric tensor perturbation (nondimensional) , connected with GW, we 

consider only the scalar, having in mind the averaged value of this perturbation. 
Taking into account h = u)h = 27rc/i/A, A ~ 10 Mpc is the wavelength of GW of 
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the order of the size of the cluster. From the comparison of (|67|) , ( 68 ) we obtain 
the averaged metric perturbation due to very long GW in the form 



/.^^f^V^'^e-io-" (69) 

C' \ TT J 

for the values of A and Sgw, mentioned above. Insidently we may expect 10 
times larger amplitude of the GW, than the averaged over the volume value. We 
may estimate the angle of the lensing using its linear dependence on the potential 
difference, and the value of deviation 1.75 arc sec. for the Sun (Zeldovich and 
Novikov, 1971) which gravitational potential is equal to hsun = GM/Rc^ w 
2 • 10"^. Finally we get for the expected lensing influence of the very long GW 

Ogw = Osun^ - 3 • lO-*0,„„ w 5 • lO-^arcsec. (70) 

Such angular resolution is available by the very long base interferometry and 
the interferometric optical telescopes under construction could give much better 
angular resolution. 



9 Conclusions 

The simplified model of the collapse and subsequent relaxation of the rotat- 
ing compressed shperoid may be related to the processes in the central parts 
ob the dark matter objects. The formation of the transient prolate dark mat- 
ter spheroid follows from the calculations for slowly rotating (or non-rotating) 
objects. 

For realistic values of the relaxation time we have obtained rather slow damp- 
ing of the oscillations, indicating that after ~ 10 oscillations the amplitude may 
remain on the level of 1/10 of the Hubble value of the velosity for a correspond- 
ing mass scale. Such oscillations in the dark matter objects of the largest scale 
may be preserved until the present time. 

Interaction of the cosmic microvawe background radiation (CMB) with the 
dark matter objects on the stage of their collapse and subsequent oscillation 
phases may influence the visible characteristics of the observed picture of the 
fluctuations of CMB. 

The weak but very long gravitationl waves (GW) emitted mainly on the 
stages of the collapse and pancake formation form a long wave GW background, 
which also influence CMB fluctuations (Doroshkevich et al, 1977). 

Besides, the existence of such a long GW in the space between the source and 
the observer may be registered as an action of the gravitational lense. Absence 
of the appropriate lensing objects in the direction of any observed gravitational 
lense object may indicate to the existence of dark matter objects without the 
barion matter presence, or may be the indication to the presence of the very 
long GW which is moving in the space between us and the lensed object The 
expected angles of such lensing (^ lO^'^) arc sec. may be available in the near 
future. 
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Figure 1: Time dependence of semiaxies a{t) and c{t), and sinus of the rotational 
phase: phase{t) for initial values of the variant 1, Table 1. 

The present approach describing the compressed spheroid by the two degrees 
of freedom may be evidently generalized for a description of the 3-axis spheroid 
by the ordinary equations with 3 degrees of freedom. In particular, it would 
be possible to investigate the influence of the damping processes on the bound- 
aries of the stability of the Maclaurin spheroid for the transformation into the 
Jacobi 3-axis ellipsoid (secular and dynamical) at large angular momentums. 
The model could be generalized also for the case of anizotropic relaxation and 
creation of anizotropic pressure. 
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Figure 3: Same as in Fig.l, for initial values of the variant 2, Table 1. 
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Figure 5: Same as in Fig.l, for initial values of the variant 3, Table 1. 
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Figure 15: Same as in Fig.l, for initial values of the variant 8, Table 1. 
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Figure 16: Same as in Fig. 2, for initial values of the variant 8, Table 1. 
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